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This letter describes a scalar curvature invariant for general relativity with a certain, distinctive 
feature. While many such invariants exist, this one vanishes in regions of space-time which can be 
said unambiguously to contain no gravitational radiation. In more general regions which incontro- 
vertibly support non-trivial radiation fields, it can be used to extract local, coordinate-independent 
information partially characterizing that radiation. While a clear, physical interpretation is possible 
only in such radiation zones, a simple algorithm can be given to extend the definition smoothly to 
generic regions of space-time. 
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The prediction of the observable properties of gravita- 
tional radiation is an important application of the the- 
ory of general relativity. There has been a recent flurry 
of activity in gravitational wave research as terrestrial 
detectors have been built to provide the first direct de- 
tection of the phenomenon. These experiments are so 
delicate, however, that detailed theoretical predictions 
are required when generating filters to identify physical 
radiation in the data. To provide these filters, researchers 
have turned to computers. An impressive array of tech- 
niques have been developed which, in the near future, 
are likely to allow dynamical simulations of astrophys- 
ically relevant systems according to the full, non-linear 
Einstein equations. Once the technical obstacles to such 
simulations have been overcome, however, tools will be 
needed to interpret their results. Here, we propose a 
tool specifically designed to extract gauge-invariant in- 
formation about gravitational radiation from space-time 
geometry. Our proposal centers on a scalar quantity con- 
structed from the space-time curvature which, in regions 
where gravitational radiation is unambiguously defined, 
manifestly contains information only about the radiation 
field. We refer to such quantities as radiation scalars. 
While the physical interpretation of our radiation scalar 
is clear only in special regions of space-time, its definition 
(generically) extends smoothly throughout. We present a 
simple algorithm to calculate it everywhere in space-time 
at only a modest computational expense. 

The notion of gravitational radiation in general rela- 
tivity does not have universal meaning. Rather, a proper 
definition can be given only in regions of space-time 
whose geometry is characterized by two distinct length 
scales: one describing an "average" radius of curvature 
and a second, much shorter scale corresponding to the 
wavelength of gravitational waves. Theoretical tools have 
been developed which allow a detailed analysis of radia- 
tion fields in such radiation zones. The best-known ex- 
amples include the standard treatment of weak waves on 
flat space-time, asymptotic techniques at null infinity jlj, 
and the Teukolsky formalism for weak radiation in the 
exterior of an equilibrium black hole M . All of these ap- 



proaches separate the physical metric into a background 
piece and a perturbation thereof describing radiation. 
They are therefore inappropriate in regions of space-time 
supporting stronger gravitational fields like those typi- 
cally modeled in numerical relativity. 

Physically, one expects a typical space-time will admit 
radiation zones both at large distances from any radi- 
ating sources, and also closer to them at late times as 
they approach a quiescent final state. However, there are 
difficulties realizing this picture in numerical relativity. 
Limitations imposed by computer technology make it dif- 
ficult to simulate evolution at large enough distance or to 
late enough times to deploy perturbative techniques un- 
ambiguously. Even within a radiation zone, it is difficult 
to separate a single, physical metric which has been cal- 
culated numerically into background and perturbation. 
There may be many flat (or Kerr) metrics on a finite re- 
gion of space-time which approximate the physical metric 
to a given accuracy at every point. Employing a radiation 
scalar can avoid these difficulties because it (a) is calcula- 
ble everywhere in space-time and (b) does not require any 
background metric be chosen. On the other hand, how- 
ever, the scalar we describe below does not offer as com- 
plete a description of gravitational waves in a radiation 
zone as the perturbative approaches mentioned above. 
In particular, it vanishes even for certain non-trivial ra- 
diation fields. However, such fields are non-generic and 
unlikely to be found by numerical simulations not explic- 
itly designed to do so. Our goal is not a formalism suited 
to all possible cases, but rather an approach which will 
be effective in generic situations one might reasonably 
expect to encounter. 

The use of scalar curvature invariants in numerical 
relativity is quite natural. Unlike the components of a 
higher-rank tensor, the value of a scalar function at any 
point of space-time is independent of any coordinates 
used to evaluate it. In a radiation zone, for example, 
a scalar can easily distinguish physical radiation from 
oscillations in the field variables induced by a peculiar 
choice of coordinates. 

Curvature invariants have been studied extensively in 
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the literature |3|, [4[ , particularly in connection to the alge- 
braic (Petrov) classification of space-times || . The best- 
known, and algebraically simplest, examples are the two 
complex quantities 
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where *C a b cd ■— \^ab mn C mn cd denotes the space-time 
Hodge dual of the Weyl tensor. In fact, these two scalars 
are not only simple, but also the only independent, alge- 
braic curvature scalars one can construct. This is because 
any curvature invariant should be expressible as a func- 
tion of a set of six fundamental (real) invariants known 
as the eigenvalues of the Weyl tensor. Since the traces 
of both the Weyl tensor and its dual vanish, a simple 
counting argument shows / and J alone suffice to recon- 
struct those six eigenvalues. Thus, any scalar curvature 
invariant should be expressible as some function of I and 
J. The operative question then is which such functions 
define radiation scalars. 

Neither quantity defined in Eq. (|l|) is itself a radia- 
tion scalar. The stationary Kerr black hole, for example, 
is known to contain no gravitational radiation anywhere 
in space-time. However, neither / nor J vanishes in the 
Kerr space-time. Intuitively, both contain contributions 
not just from the radiative part of the gravitational field, 
but also from the Coulombic part. We must therefore 
seek a non-trivial function of both / and J in which their 
Coulombic contributions cancel. One can adopt a sys- 
tematic approach to this search. 

The Newman-Penrose formalism || encodes the ten 
independent components of the Weyl tensor at each point 
of space-time in the five complex scalars 
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These curvature components $„ depend on one's choice 
of a null tetrad (£ a , n a , m a , rh a ) on space-time. All four 
vectors in such a tetrad are null, with the first pair real 
and the second both complex and conjugate to one an- 
other, and their only non-vanishing inner products are 
l a n a — — 1 and m a rh a = 1 . The Newman-Penrose ap- 
proach is particularly useful in the search for radiation 
scalars because it already expresses space-time curvature 
in the coordinate- independent scalars *& n . However, this 
independence comes at a price. In general, there is no 
preferred null tetrad on space-time, and clearly differ- 
ent choices in Eq. (|^) will lead to different values of 
There are then two ways to extract invariant, physical re- 
sults. One can either seek combinations of the \t n which 
take the same value in every frame (such as / and J), 
or find a means of "gauge-fixing" by imposing additional 
conditions on the tetrad. We shall employ a combina- 
tion of these techniques below. We should emphasize 



that, while the Newman-Penrose formalism is a conve- 
nient tool in our search, the final result we present below 
does not depend on it. Rather, the radiation scalar we 
propose can be calculated directly from the space-time 
curvature; it is a particular function of I and J. 

In a local reference frame adapted to a given null tetrad 
(£ a , n a , m a , rh a ) , the ^ n have specific physical interpre- 
tations 0. The components \Po an d (^4 and ^3) 
describe transverse and longitudinal gravitational waves, 
respectively, propagating in the null direction £ a (n a ), 
while ^2 represents the Coulombic part of the field. 
These interpretations are supported by results in the 
weak-field limits associated with radiation zones. More- 
over, both near null infinity and in the Teukolsky formal- 
ism, physical radiation is associated with the transverse 
components ^0 an( i ^4 while the longitudinal compo- 
nents and ^3 are understood to be "pure gauge." In 
the Teukolsky formalism, for example, there is a natural 
null tetrad. It is the Kinnersley frame , whose real null 
vectors lie in the two repeated principal null directions 
of the background Kerr geometry. The null tetrad in the 
Teukolsky formalism differs from the Kinnersley frame 
only at first order in perturbation theory. However, the 
exact perturbation in the tetrad can always be chosen 
such that 
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Motivated by these considerations, we seek to define a 
radiation scalar which is a function solely of ^0 an d x f , 4- 
However, the apparent simplicity of this goal is mislead- 
ing. It only acquires meaning if one can find a null tetrad 
everywhere on space-time which passes over smoothly to 
the correct limits at late times and large distances. 

While the condition Eq. (||) has been introduced in 
the context of perturbation theory on Kerr space-time, 
there is no obvious problem imposing it more generally. 
We shall call a null tetrad satisfying these conditions a 
transverse frame. It is not immediately clear whether 
transverse frames always exist and, if they do, whether 
they are unique. However, it is possible to answer both 
of these questions locally at any point of any space-time. 

When they do exist, transverse frames are never 
unique. This is because there are two types of null tetrad 
transformation which preserve the conditions of Eq. (|^). 
These include a continuous family of spin-boost transfor- 
mations S(c) which rescale both real null tetrad vectors 
and rotate the complex ones by a phase, as well as a 
discrete exchange transformation E which interchanges 
the real and complex vectors in pairs. They affect the 
curvature components according to 



S(c) : i-> c 2 - n and E : 
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These ambiguities in a transverse frame always exist, 
however, and may not exhaust one's freedom in choos- 
ing it. It is therefore natural to ask how many distinct 
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transverse frames up to spin-boost and exchange exist at 
a given point of space-time. This question has a remark- 
ably simple answer. Each transverse frame is associated 
with a non-null, self-dual eigenvector of the Weyl ten- 
sor. This association is not one-to-one, but rather two 
transverse frames give the same eigenvector if and only 
if they can be transformed into one another using only 
spin-boost and exchange operations. While the details of 
the association are unimportant here, it gives the crucial 
result that the number of distinct transverse frames is 
completely determined by the algebraic (Petrov) class of 
the Weyl tensor. In the various cases, one finds 

Type I There are exactly three transverse frames. 

Type I-D There are infinitely many transverse frames, 
falling in two sets. One set contains only the Kin- 
nersley frame, while the other contains a continuum 
of other transverse frames. 

Type II There is only one transverse frame. 

Type II-N There are a continuum of transverse frames. 

Type III There are no transverse frames. 

Geometries of types II-N and III describe pure radiation 
fields, and are rarely encountered in numerical relativity 
unless explicitly sought. The absence or surplus of trans- 
verse frames in these cases is therefore of little practical 
consequence. Moreover, a type I-D geometry's Kinners- 
ley frame is singularly important physically. Thus, in all 
cases of primary physical interest, transverse frames do 
exist. The relevant ones are unique up to spin-boost and 
exchange, modulo an additional finite ambiguity in the 
type I case which we resolve below. 

In the absence of additional structures, such as a 
Killing field, one cannot proceed to eliminate the resid- 
ual spin-boost and exchange ambiguities in a transverse 
frame. We therefore adopt a different approach. Accord- 
ing to the transformation laws, Eq. (|]), the quantities 

€ := *o *I ( 5 ) 

take the same value for all transverse frames related to 
one another by residual transformations. (The notation 
Vf^ indicates evaluation in a transverse frame.) In a type 
I geometry, there will be three distinct, non-zero values of 
£, while for type II the only value is £ = 0. There will be 
two distinct values in a type I-D geometry. One, £ = 0, is 
associated with the Kinnersley frame, and the other with 
the continuum of other transverse frames. Motivated by 
all these facts, we seek a function on space-time which 
equals one of the values of £ at every point. If this func- 
tion always approaches zero in the limits of late times 
and large distances, it will define a radiation scalar. 

At first, the prospect of evaluating £ numerically seems 
daunting. Indeed, its definition involves one or more alge- 
braically preferred null tetrads which perhaps cannot eas- 




FIG. 1: The Riemann surface Zs is a triple-cover of the Rie- 
mann sphere. The values of Z(S) on its three sheets corre- 
spond to the three values of £ . The labels in the figure indicate 
contiguous regions of Zs pictured on different sheets. 

ily be found numerically at each space-time point. How- 
ever, the close connection between transverse frames and 
the algebraic structure of the Weyl tensor suggests an al- 
ternative. One can construct a tensor S a h cd , a quadratic 
polynomial in C a b cd , whose eigenvalues are precisely £. 
The characteristic equation for this tensor gives 

P := e - + ^ - ^ (/ 3 - 27 J 2 ) = 0. (6) 

The three roots of P are the values of £. This result shows 
explicitly £ can be evaluated directly without invoking 
the Newman-Penrose formalism. 

The roots of Eq. @ take a particularly simple form 
when expressed in terms of the Baker-Campanelli spe- 
ciality index |5) S := 27 J 2 // 3 . This quantity equals 
unity if and only if space-time is of types II or I-D (it is 
ill-defined for types II-N and III, where I = = J). The 
classical Cardano formula yields 

C = IZ{S) := ~ (2- [W(S)} 1/3 - [W(S)}- 1/3 ) , (7) 

where W(S) := 2S - 1 + 2\/S 2 - S. The three values of 
the cube root yield the three solutions of Eq. (pi) . 

The function Z(S) is not single-valued, and one must 
therefore pick a branch of it to evaluate £ numerically. A 
careful analysis shows the Riemann surface Zs 1 depicted 
in Fig. [l], for Z(S) consists of three sheets with branch 
points at 0, 1, and 00. There is a unique sheet Z$ with 
no branch point at S = 1. The associated branch Zq(S) 
is unique among the three in that Zq(1) = 0, whence only 
this branch gives the correct limit, £ ~ 0, in the radiation 
zones at late times or large distances. Thus, a function 
which equals £0 : = I Zq{S) in these limits will define a 
radiation scalar. One can evaluate £0 numerically using 
branches of the roots in Eq. ([?]) with \/l = VI = 1 , and 
branch cuts for both along the negative real axis. With 
these choices, the branch cut for Zq(S) itself also lies 
along the negative real axis. 

A potential problem arises as one attempts to extend 
£0 into regions of space-time supporting stronger gravita- 
tional fields. Namely, at certain points along a curve con- 
necting a well-understood radiation zone to the strong- 
field interior, the value of S may lie on the branch line 
for Zq(S). At such points, £0 as defined above will be 
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discontinuous. However, although a given branch of a 
multi-valued complex function is necessarily discontinu- 
ous at a branch line, it is always continuous there with 
one of the other branches. In our case, this means we 
can define a function along the curve in question which 
equals £o for points with values of S just to one side of the 
branch line, and equals one of the other roots, say £+, for 
nearby points with S just to the other side. This function 
is guaranteed to be smooth on space-time, and is simple 
enough to evaluate numerically. Begin by calculating all 
three roots of Eq. (^j) at each point of space-time. At 
late enough times and large enough distances, pick the 
root £o- Then, extend the function globally by moving 
to points with successively stronger fields, always picking 
the appropriate root so as to keep it smooth. This algo- 
rithm, for a generic space-time, will result in a smooth 
function whose large-distance and late-time limits will 
equal £o- It is a radiation scalar. 

The speciality index is a convenient measure of 
whether a simulation has reached far enough into the 
asymptotic regions for the root £o to be clearly rele- 
vant. Namely, if one finds in numerical data a fairly 
large boundary region where \S — 1| < 1, it becomes rea- 
sonable to assume curves extending outward to infinity 
will not support values of S lying on the branch line for 
Zq (S) . Practically speaking, this gives a quantitative in- 
dication of the whether the above construction can be 
meaningfully applied to a given simulation. 

One might worry the radiation scalar we have de- 
scribed is not single- valued. Indeed, it is defined by 
extending its value smoothly along curves linking well- 
understood radiation zones to other, strong-field regions 
of space-time. In principle, the value at a point in the 
strong-field region could depend on the curve used to 
reach it. However, for the vast majority of space-times, 
this worry is unfounded. Specifically, if the algebraic class 
of space-time is nowhere exactly type II or higher (II-N 
or III; there is no problem in type I-D), the radiation 
scalar is necessarily single-valued. It may be so in other 
situations as well, but a proof does not exist. 

A conspicuous feature of the definition of £ in Eq. (|J) is 
its dependence on both \&o and ^4. Near infinity, these 
curvature components describe the ingoing and outgo- 
ing parts, respectively, of the radiation field. There are 
therefore non-trivial (pure ingoing or outgoing) radiation 
fields for which £ will vanish. However, such configura- 
tions are not generic. Real radiation fields experience 
some back-scattering, even at large distances, and only 
for very carefully designed states will ^0 or ^4 vanish ex- 
actly. In fact, the ingoing and outgoing parts of the radia- 
tion field are coupled through the Teukolsky-Starobinsky 
identities |l(], [ll| ; complete knowledge of either one de- 
termines the other. Thus, for example, insisting the in- 
going part vanishes exactly imposes severe restrictions 
also on the outgoing part. As a result, pure ingoing or 
outgoing fields are more difficult to construct, and rarer, 



than one might first expect. 

As a sample application of these techniques, consider a 
full, non-linear simulation of a binary black hole system. 
After coalescence, as the single remaining hole settles 
down to a quiescent final state, the relevant transverse 
frame will approximate the Kinnersley frame of some (un- 
determined) background Kerr geometry. The factors, *q 
and defining £0 will therefore satisfy the Teukolsky 
equation. Their late-time behavior is known from pertur- 
bation theory: after quasi-normal ringing, both factors 
decay to zero as a power of time. Moreover, the (com- 
plex) quasi-normal frequencies are specific functions of 
the mass and angular momentum of the final black hole. 
One can therefore extract these physical parameters from 
the late-time behavior of the radiation scalar without first 
invoking perturbation theory techniques. This is impor- 
tant conceptually since the ideal background for a given 
simulation would have the physical parameters so deter- 
mined. They would not otherwise be known when mak- 
ing the transition to perturbation theory. 

While it may be possible to construct radiation scalars 
other than the one described here, they may not be easy 
to find. A natural candidate, g := S — 1, has many fea- 
tures similar to those of £0 , but does not have the correct 
behavior near null infinity. To leading order in pertur- 
bation theory, this quantity is given by ipo V'Z/W'J ) 2 - 
While g certainly vanishes if space-time is exactly Kerr, 
the denominator here generally falls off at the same rate 
as the numerator. Although the asymptotic value will be 
small, quadratic in the perturbation parameter, it will 
not fall off to zero. It therefore does not extract informa- 
tion solely about radiation near null infinity. 
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